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Abstract 

We calculate the atomic (spin) Wigner function for the single mode Dicke model in the regime 

>' 

. of large number of two-level atoms. The dynamics of this quasi-probability function on the Bloch 

(N 

sphere allows us to visualize the consequences of the entanglement process between the boson and 

' the spin subsystems. Such investigation shows a distinct localization behavior of the spin state 

with respect to the polar and azimuthal Bloch sphere angles. A complete breakdown of reflection 

Q-i. symmetry in the azimuthal angle is shown in the non-integrable case, even at short evolution times, 

-i— > ' 

Also, in the classically chaotic situation, the appearance of sub-planck structures in the Wigner 

C3 ' 
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X 



function is shown, and its evolution analyzed. 
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I. INTRODUCTION 



The Wigner function (WF) [1] has been an important tool in physics, particularly to ex- 
plore the interface between classical and quantum world [2]. Furthermore, it is an useful 
ingredient in the theory of partial loss of coherence [3]. Recently, advances in the state 
reconstruction problems [4] renewed interest on the WF, particularly after the measure- 
ment of negative quasi-probabilities [5, 6] which are generally considered as signatures of 
non-classical states [7]. Also, quasi-probability functions are relevant for the study of en- 
tanglement process of continuous variables, as in the beam splitters [8], in order to analyze 
separability of two-mode Gaussian states [9, 10, 11], and in studies of teleportation of non- 
classical states like in Ref. [12]. Since the proposal of manifestation of chaos in the process 
of decoherence by Zurek and Paz [13], the WF has also been an important tool for studying 
this issue [14, 15]. 

Here, our focus is the use of the atomic Wigner function to study the ideal atom-field 
entanglement process. Recent experimental progress on arrays of quantum dots and Joseph- 
son junctions has raised the possibility of super-radiance, a phenomenon already known to 
arise in a system with N two-level atoms coupled to a field mode described by the Dicke 
model (DM) [16]. Some of the solid-state systems relevant to quantum information that can 
be mapped into DM are: 'phonon cavity quantum dynamics' [17, 18]; those with Joseph- 
son Junctions and quantum dots, with the possibility of multiqubit entanglement Dicke-like 
model [19]; and the proposal of 'circuit QED' [20], which has been constructed recently for 
one qubit Jaynes-Cummings case [21]. Finally, the atom-molecule coexistence model near 
a Feshbach resonance of cold fermions in both strong [22] and weak [23] coupling has been 
proposed as a DM. 

We are interested in exploring the atomic (spin) Wigner function dynamics of the gen- 
eralized N- Jaynes-Cummings model (A-JCM) [24] where both Jaynes-Cummings (JC) and 
anti- Jaynes-Cummings (AJC) interaction can be present. These interactions between spin 
and boson systems has been shown to appear naturally in trapped ions [25, 26], or in cavity 
QED by means of strong classical driving field [27]. Classically the model becomes chaotic 
[28, 29, 30], and some of its manifestation in the entanglement process has been shown in 
previous publications [31, 32, 33, 34, 35]. To our knowledge, it is the first time the dy- 
namical evolution of the A-qubits atomic Wigner function is followed in detail during the 
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entanglement process since it has been introduced [36, 37, 38]. 

Remarkable features in the AWF can be seen in such wave-packet dynamics, showing 
dynamical differences between localization of the angular variables 0, 9 on Bloch Sphere. We 
also explore the dynamics of the negative-valued parts of the AWF, usually considered as a 
hallmark of non-classicality and interference effects. Another notable aspect is the distinct 
time evolution of the AWF, depending on the initial positions chosen for the centers of the 
wave packets. Such sensitivity to initial conditions have been already noted in references 
[31, 32], at the level of the integrated quantities like the subsystem linear entropy or mean 
values, but the details of the behavior of quasi-probability functions was lacking. 

This paper is organized as follows: in Sec. II we present the classical analog of the Dicke 
model to define and show some quantities used in the following sections. Also, we review 
the notion of sensitivity of the entanglement process on the initial position of the coherent 
packet, in both integrable and chaotic cases, already found in previous works. Sec. Ill is 
reserved to present our results for the A-qubit atomic Wigner function for some selected 
initial conditions and analyze its time evolution under the point of view of the entanglement 
process. In Sec. IV we summarize our results. 

II. THE MODEL AND ENTANGLEMENT DYNAMICS WITHIN THE LARGE- N 
APPROXIMATION 

We consider a generalized version of A-atom (qubit) Dicke model, with variable coefficient 
for the rotating and counter-rotating wave term. Moreover, since our aim here is to pursue 
the effect of the classically chaos generating term on the entanglement between the A-atom 
and the field, we will ignore the interaction between the atoms (qubits) and treat the atomic 
system as a large spin (A = 2 J). 



Here, ujq and uo a are frequencies associated with free Hamiltonians for field and atoms respec- 
tively. G, G' are coupling constants associated with atom-field interaction within the dipole 
approximation. The usual Rotating Wave Approximation (RWA) is recovered by setting 
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G' — 0. The field observable are described by means of the creation and annihilation opera- 
tors a and a\ whereas J z , J± are pseudo-spin operators associated to an atomic observable. 
This model is used to describe both, cavity QED experiments [27, 39] (with G = G", but 
usually to an excellent approximation one can set G' = 0) and trapped-ion systems. In the 
last system, interactions with different couplings G ^ G' can be generated [25, 26]. 

For the purpose of our study where both systems are initially pure and separable in 
quasi-classical states, the appropriate initial state \wv) is a product of the field and atomic 
coherent states defined as [36, 40, 41]: 



M 2 -t 



\u) = D(v) |0) = e~^e ua +u a |0) 
\w) = (l + \w\ 2 )~ J e wJ+ \J,-J). (2) 

Here, J = N/2 and the variables w and v can be written as a function of the classical 
variables in the corresponding phase spaces, (?/,£>/) for the field, and (q a ,p a ) f° r the atomic 
degree of freedom 

_ Pa + Mia 

y/4 J - (pl + ql) 
v = (Pf + Wf) (3) 

A corresponding classical Hamiltonian can been obtained by a standard procedure [42] , using 
the above defined coherent states (wv\ H \wv) [30] 

U (q a ,Pa, qf,Pf) = ^{ P } + qj) +^-{pl + Ql) (4) 



4 J - [pj + q 2 f/ 
~u a J + 1 —j (G + p a p f + G_q a q f ) , 

with G± — G ±G'. The classical dynamics associated with this Hamiltonian were explored 
before [30], and shown that:(jj integrable situations are recovered when either G or G' is zero; 
(ii) the most chaotic dynamics is associated with the condition G = G' ~ O(ivo,iv a ), as we 
increase the coupling constant. Our aim is to investigate the time evolution of the initially 
quasi-classical wave-packet along with the occurrence of entanglement between the iV-atom 
and the field systems. Particularly, we are looking for the possible differences in the reduced 
wave-packet dynamics when we compare integrable and chaos generating interactions. The 
connection with the classical dynamics is established by choosing coherent states as initial 



states, centered at the corresponding points of the phase space. Then, we let the system 
evolve by means of the Hamiltonian (1) and explore the entanglement dynamics solving 
numerically the Schrodinger equation. 

In order to know where to put the initial atomic and field wave-packets, the first step is 
study the structure of classical phase space. This can be done by using the Poincare section. 
In Fig. 1, we plot the projection of Poincare section in the atomic plane (q a ,p a ) for qf = 0. 
and pf > 0. The surfaces of section correspond to both, the right one to the integrable 
case with (G = 0.5, G' = 0) and the left one to the soft chaos (G = 0.5, G' = 0.2). Here, 
and along the present work, the total energy is fixed at E = 2J = 21.0 with J = 10.5. 
The coupling values corresponds to the non-super-radiant phase {G + < 1) [43]. The limit 
of atomic phase space is indicated by a border at radius equal to \/AJ . Integrable section 
shows a separatrix of motion along the line p a = 0.0 and concentric tori around each of the 
two stable periodic orbits. A pro-eminent feature of the non-integrable surface of section is 
a large stability island for p a > 0. 

The symbols in Fig. 1 show the chosen centers for the atomic coherent states. We choose 
two specific initial conditions (i.e.) for each case. In the integrable case, the first one (II) 
is on an internal tori belonging to the region p a > 0, marked by a triangle; and the second 
(12), marked by a circle, located near the border of the atomic phase space. For soft chaos 
situation, the first i.e. (Nl) where chosen on a point inside the largest stability island (circle). 
The second condition (N2) is located in the chaotic sea (triangle). Specific values of q a , p a 
and the mean value of J z operator for the associated atomic coherent state are listed in 
Table II. 
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(Jz)/J 


11 


0.0 


0.55 


-0.43 


12 


0.1 


0.95 


0.83 


Nl 


0.0 


0.54 


-0.47 


N2 


0.0 


-0.28 


-0.84 



TABLE I: q a , p a and (J z ) values for the chosen initial conditions. 

As a second step, one should calculate the time evolution of the atomic linear entropy 
(ALE) which, in this globally pure bipartite system, can be used as a measure of entangle- 
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FIG. 1: Poincare section for the atomic degree of freedom (with qj = 0.0 and pf > 0) in the 
resonant case uj a = and energy E = 21., with N = 2 J = 21. Left: integrable case with G = 0.5 
and C = 0.0. Right: non-integrable case with G = 0.5 and G" = 0.2. 

ment. This calculation is done in three stages: First we diagonalize the Hamiltonian (1) 
numerically and use the eigenvalues and eigenstates in order to find the temporal total 
density operator, pit). Then, we calculate the reduced density operator of the atomic sub- 
system, p a = Tif [p (t)], and finally obtain 5 a (t) — 1 — Tr a p a 2 (t) . Evolution in time of the 
ALE for i.c.'s II and 12 are plotted in Figure 2(a). For both i.c.'s, the subsystem entropy 
increases in the mean as time goes on, until a plateau is reached. The details of the curve 
such as the specific values of the ALE in the plateau, the particular oscillatory behavior 
and the entanglement rate depend on each initial condition. In particular, we observe that 
those atomic i.c.'s with the classical dynamics restricted to a well delimited region in phase 
space on the tori region are more resistant to entangle with field. This relation between fast 
entanglement process and less localized classical dynamics was pointed out in a previous 
work [32]. 

The ALE for the soft chaos situation are plotted in Fig. 2(b). In this plot, it is clear that 
larger entanglement rate is also associated to the chaotic i.c.'s. It is interesting to compare 
our results for II and Nl, solid lines in Fig. 2(a,b). For II, ALE shows more regularity in 
the oscillations and reaches the corresponding plateau around t ~ 30, whereas in ALE for 
Nl the oscillations are less regular and takes a longer time to reach the plateau, at t m 70. 
Also, ALE for Nl condition keeps a small oscillation around its mean value (due to the 
non-RWA term) and one can also see a certain large period modulation in contrast to II. 
The relation between the behavior of the classical trajectory and maxima and minima of S a 
was studied in Ref. [33], and we will not discuss here. Instead, the atomic Wigner Function 
(AWF) is used to visualize the behavior of global state as projected in atomic phase space. 
In the next section, we will illustrate the following aspects: (i) that the AWF provide us 
the information that those i.c.'s that are more "protected" against entanglement process, 
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FIG. 2: Atomic Linear Entropy of the iV-JCM associated with the four initial conditions centered 
at the positions shown in Fig. 1 (with the same parameter values) and listed in Table II. (a) 
Integrable case: i.e. II (Solid line) and 12 (dotted line); (b) Soft chaos: i.c Nl (Solid line) and N2 
(dotted line). 

have a strong localization on the Bloch sphere; (ii) we also show how the destruction of tori 
due to chaos in the classical dynamics goes along with the derealization of the quantum 
wave-packet during its temporal evolution. 

III. DYNAMICS OF ATOMIC WIGNER FUNCTION. 

The phase space quasi-probability distributions of electromagnetic field and atom have 
been discussed by several authors [1, 36, 37, 44, 45]. Here, we adopt the definition of Wigner 
function in terms of arbitrary angular momentum basis, as introduced by Agarwal [37]. This 
function is defined as 



This is the characteristic function associated with atomic density operator p a (t). Here, Tkq 
is the multipole operator acting in the angular momentum space [46] 




k=qq=-k 



(5) 



where Qk,q is given by 



Qk,q (*) = Tr p a (t) T K q 



(6) 




x |J,M) (J,M-Q\. 



(7) 
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In Eq.(5), the usual Wigner 3 J symbol has been used, and Y K q (9, 0) indicates the spherical 
harmonics defined over the Bloch sphere. There, 9 is the polar angle and is the azimuthal 
angle. The distribution of any angular momentum state can be studied using the AWF. As 
shown in Dowling et al. [38], it is possible to estimate the indeterminacy in the measure of 
the atomic observable J z , J x J y in the state through its variances . In fact, if the state has 
a large probability associated to a well-defined eigenvalue of J z , its atomic Wigner function 
shows a strong localization in polar angle 9. In a similar way, indeterminacy associated to 
the measure of J x and J y means ignorance on the azimuthal angle 0. 

In order to obtain the AWF, we use our previous results of atomic density matrix operator. 
Because the basis used was the Dicke states | J, J z ), we can calculate the action of Trq on 
each atomic density matrix elements, obtaining W (8,<f),t). Also, we always set the = 
value exactly at the center of each initial atomic coherent packet. That means, if the wave 
packet is not on the X-axis, we rotate the XY plane by a certain angle 0(t = 0) = O in 
such a way that the direction defined by the vector (sin0 o , — cos0o, 0) coincides with the 
rotated X-axis. 

A. Integrable case with G = 0.5 and G' = 0. 

For the integrable case, we show the snapshots of the temporal evolution of contour lines 
of AWF in Fig. 3 for the i.e. II. The initial coherent state, Fig. 3(a), has its maximum 
value (3.5) at 9 = 0.647T, which corresponds to (J z ) « — 0.43J. At the time when the first 
maxima of ALE is reached, the atomic state has a more delocalized distribution, shown in 
Fig. 3(b), with two negative valued regions (in black). A formation of three overlapping 
positive peaks starts, with maxima (almost) at the equator of Bloch sphere (9 = 7r/2). At 
this time, AWF has a maximum value lower than the one at the initial time (~ 1.7). The 
appearance of a negative part in the AWF with value —0.2 (10% of maximum value) 
indicates the non-classical character of this state. 

The derealization of the state in the azimuthal angle is associated with the increase of 
the ALE. This assertion can be confirmed by checking the evolution for consecutive maxima 
and minima. The forms of AWF are shown at times corresponding to the first minimum, 
Fig. 3(c); second maximum, Fig. 3(d), and second minimum, Fig. 3(e). Comparing them, it 
is clear that the AWF is more localized in the variable at times when ALE has minima. 
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Although in Fig. 3(d), the three peaks have coalesced into one, the packet sweeps a larger 
interval over values than those found at the first minimum time scale. We can also 
observe that AWF maximum value oscillates around 9 xs 1.8. A negative valued portion is 
still present but it became significantly smaller (less than 1%) than in Fig 3(b), so it is not 
possible to see in Fig 3(c,d), but reappears in Fig. 3(e). From this sequence, it is clear that 
atomic state looses both, azimuthal and polar localization, associated with the increase in 
the atomic linear entropy. However, at those times when the ALE plateau is reached, we 
see how the AWF (plotted in dotted lines) still has localization in 9 (near the equator in the 
Bloch sphere), having non-zero values only in the interval 1.6 < 9 < 2.5. This shows that for 
the particular initial condition considered on an internal torus (distant from the separatrix 
motion and the border of the phase space), an increase in the ALE is associated with the 
increasing derealization with respect to the azimuthal angle much more than to the polar 
variable. 

Now we shall compare previous case with the time evolution of the AWF for condition 
12 (9 ~ 0) which is close to the largest value of (J 2 ), shown in Fig. 4. At the initial time, 
Fig. 4(a), the AWF is well localized in both angular variables and, as the two subsystems 
interact, we can see that the AWF begins to lose localization mostly in the polar variable 
(non-zero values in the interval 1.0 < 9 < 2.7), at the time scale of the first maximum in 
ALE. Physically, the spreading of the AWF and the increase in the atomic linear entropy, 
tell us that for this i.e. (12), the increasing entanglement of the iV-atoms with the field 
is mostly associated with an increase of the participating states of atomic levels in the J z 
spectrum (shown by the spreading of the AWF in 9 variable). The value of the ALE at 
the first maximum for 12 (5 a ~ 0.5) is greater than for II (5 a lower than 0.2), where the 
increasing derealization is primarily in the azimuthal variable <fi. 

At the time corresponding to the first minimum of the ALE, Fig. 4(c), we observe that 
AWF forms a three-peak structure seen in the II case. However, there is a appreciable 
spreading in the variable but not allowing the peaks to become separated, and near the 
central position in <fi = occurs also a notable derealization in 9 variable. The negative part 
is no longer localized as in the previous i.e., and its value is ~ 10% of the maximum of AWF. 
Again, we observe the connection between minima of 5 a and the angular localization: the 
AWF for the second maximum of ALE, Fig. 4(d), shows how the atomic state is completely 
delocalized in both variables at this time. The AWF is non-zero practically at all points 
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FIG. 3: Shaded contour plots of atomic Wigner function associated with the i.e. II: (a) initial 
time; (b) first maxima of the corresponding atomic linear entropy (ALE); (c) first minimum of ALE; 
(d) second maxima; (e) second minima; (f)i = 40 at the plateau region. Negative valued regions of 
AWF are drawn in black. 

on the atomic phase space. This situation is reversed at the time of the second minimum 
of ALE, with a tentative to re-gain some localization in 6 variable and a structure which 
roughly resembles Fig. 4(b). It is also interesting to see how the negative part of AWF 
reappears and it is even more pronounced than in the previously referred time (~ 20%). At 
times when the plateau is reached the AWF is totally delocalized and not even a signal of 
a main positive peak is present, which was the case we found in the plateau times of the 
internal torus case. Negative part is less than 1% of maximum value of AWF. 

From this results, it is clear that initial condition II dynamically preserves the localization 
of the AWF, especially in the polar angle. This is associated with a certain inhibition in the 
entanglement process. The difference on the value of the atomic linear entropy between the 
two initial conditions presented here is clearly related with the derealization process in both 
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FIG. 4: Shaded contour plot atomic Wigner function associated with the i.e. 12: (a) initial time; 
(b) first maxima of corresponding atomic linear entropy (ALE); (c) first minimum of ALE; (d) 
second maxima; (e) second minima; (f)i = 40 at the plateau region. Negative valued regions of 
AWF are drawn in black. 

azimuthal and polar angles. Hence, the dynamics of the internal tori is protected against 
the entropy increase, and this is related with the localization in the polar angle; whereas, 
the i.c.'s located near the separatix and the border do not have such a dynamical protection. 
Other initial coherent states with similar characteristics has qualitatively analogous behavior 
for the ALE. 

Another characteristic is the clear appearance of some sub-planck structures, namely the 
structures with considerable amplitudes with their supports in areas much smaller than h 
in phase space, similar to those discussed by Zurek [15]. On the Bloch sphere, the minimum 
action area (h is taken to be 1 here) is associated with the size of atomic coherent state 
(at t — 0) which defines a minimum-uncertainty packet. This minimum action area can be 
inferred, for instance, in Fig. 3(a) and Fig. 4(a) for N=21. Such sub-planck structures which 
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are peaks confined in areas significantly smaller than the size of the initial packet in the 
((/>, #)-plane appear for instance in Figs 4(f), 5(d,f) and 6(f). It is interesting to note that its 
appearance is indeed connected with the time where the atomic subsystem has lost its own 
coherence by entangling with the field. What is remarkable is that, for the integrable case, 
the entanglement process leaves the AWF with sub-planck structure only for i.e. 12 but not 
for II. This is, to one side, AWF counterpart of the rapid loss of coherence that occurs 
for the wave packet located near the separatrix of motion [32], but it is more than simply 
accelerating the entanglement process: the dynamical instability also generates structures 
similar to chaotic case as we shall see in the next subsection. 

B. Non-integrable case: G = 0.5 and G' = 0.2 

Now, we present in Figures 5-6 our results for the time evolution of the AWF corre- 
sponding to the non-integrable case for the conditions Nl (regular region) and N2 (chaotic 
region). Some similarities between integrable and non-integrable cases can be noticed: first, 
the connection between oscillatory behavior and a delocalization-localization of AWF are 
still present even in the chaotic i.e. (N2). This can be seen, for example, in Figs. 6(c, e), as 
compared with Fig. 6(d). There, AWF seems to suffer a "recoil" to a restricted area in (0, 9) 
plane at times which correspond to a minimum in the ALE (corresponding to the behavior of 
Figs. 4(c, e) compared with Fig. 4(d) of the integrable case). Second feature is related with 
the similar "tori protection" that was found in the integrable case. Comparing the forms 
of AWF in the plateau region, Fig. 5(f) and Fig. 6(f), a more localized AWF (particularly 
for the positive-valued part) is evident for the first condition (Nl) inside the large stability 
island than the second one (N2) in the chaotic region. Notice that, we also obtain a certain 
difference for the ALE plateau values in Fig. 2(b). 

The most interesting aspects are the dynamical differences between the two cases. Notice 
that, in the integrable case, AWF has an azimuthal mirror symmetry: <fi — > — 0, the < 
region being a mirror image of the <f) > region. This symmetry is not present in the non- 
integrable case. At this point, it is important to recall that this symmetry breaking is already 
observed in the classical Poincare section in atomic phase space (see Fig. 1). Since we are 
studying a situation within the large- N limit, this symmetry breaking can be associated with 
quantum chaos at the level of spectral distribution [29]. Another distinguishable feature is 
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FIG. 5: Shaded contour plot of the atomic Wigner function associated with i.e. Nl: (a) initial 
time; (b) first maxima of corresponding atomic linear entropy (ALE); (c) first minimum of ALE; 
(d) second minimum; (e) third maxima; (f) t = 70 at the plateau region. Negative valued regions 
of AWF are drawn in black. 

the behavior of the sub-planck structures in the AWF. They have appeared in Fig. 5(b), in 
spite of the tori protection, and remain for times at the plateau region. The size of this sub- 
planck structures seems to saturate after the entanglement time, confirming for the present 
model the results shown by Zurek. Also, it is notable the presence of a larger number of 
negative sub-planck packets than positive ones, although we do not have any explanation 
for this fact. 



IV. SUMMARY 

This work gives a complete analysis of the temporal behavior of the entanglement process 
in the N-JCM in the large- N wave packet dynamics. Previous results have pointed out the 
sensitivity to initial conditions of the atomic linear entropy. Here, a calculation of the atomic 
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FIG. 6: Shaded contour plot of the atomic Wigner function associated with i.e. N2: (a) initial 
time; (b) first maxima of the corresponding atomic linear entropy (ALE); (c) first minimum of 
the ALE; (d) second maxima; (e) second minima; (f) t = 70 at the plateau region. Negative valued 
regions of AWF are drawn in black. 

Wigner function allowed us to uncover additional information about the atomic subsystem, 
not visible in an integrated quantities like the entropy. This allowed us to have a better 
idea of what is happening to the reduced atomic state during the entanglement process, as 
a function of both, the type of interactions present (rotating and counter-rotating) and the 
initial position of the coherent wave-packet. 

A very conspicuous information obtained in this way, is the dynamics of the amount of 
the derealization of the AWF during the entanglement process as a function of both angular 
variables on the Bloch sphere. Also, we show that the presence of the classical tori structure 
in the phase space surrounding the center of the coherent wave packet, is an indication 
at the quantum level of a certain inhibition in the coherence loss. Thus, by breaking the 
integrability we also break this protection against derealization in the polar angle for the 
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initial coherent state centered on the internal tori. However, the regular surrounding is still 
an indication of slower loss of coherence and, apparently the larger the island of stability 
around the wave packet, stronger is this effect on the quantum wave packet. The symmetry 
breaking of the AWF for any time (t > 0) is another characteristic of the non-integrable 
case. 

The most interesting aspect is the dynamics of the sub-planck structures: it is completely 
absent in the regular initial conditions of the integrable situation, but do appear in the long- 
time (ALE plateau region) behavior for the packet placed near the separatrix of motion. In 
the non-integrable case, even the wave packets placed inside the regular island do develop 
sub-planck structures well before the plateau of the ALE is reached, and the chaotic cases 
show such structures already at the first maximum of the ALE. Such sub-planck structures 
seems to be directly associated with the destruction of the "tori protection" and faster 
entanglement due to less restricted dynamics in phase space, thus being an indicator of 
dynamical instability connected to quantum chaos. 
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